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Abstract 

The strength distributions of the giant monopole resonance (GMR) have been measured in the 
even-A Sn isotopes (A=112-l24) with inelastic scattering of 400-MeV a particles in the angular 
range 0°-8.5°. We find that the experimentally-observed GMR energies of the Sn isotopes are 
lower than the values predicted by theoretical calculations that reproduce the GMR energies in 
208 Pb and 90 Zr very well. From the GMR data, a value of K T = -550 ± 100 MeV is obtained for 
the asymmetry-term in the nuclear incompressibility. 

PACS numbers: 24.30. Cz; 21.65. +f; 25.55.Ci; 27.40. +z 
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Incompressibility of nuclear matter remains a focus of experimental and theoretical in- 
vestigations because of its fundamental importance in defining the equation of state (EOS) 
for nuclear matter. The latter describes a number of interesting phenomena from collec- 
tive excitations of nuclei to supernova explosions and radii of neutron stars The Giant 
Monopole Resonance (GMR) provides a direct means to experimentally determine the nu- 
clear incompressibility. 

Experimental identification of the GMR requires inelastic scattering of an isoscalar 
particle-the a particle, for example-at extremely forward angles, including 0°, where the 
cross section for exciting the GMR is maximal. Such measurements have improved con- 
siderably over the years and it is now possible to obtain inelastic spectra virtually free of 
all instrumental background directly and in coincidence with proton- and neutron-decay 
[sj]. In recent work, the GMR strength distributions have been extracted in a number of 
nuclei from a multipole-decomposition analysis (MDA) of such "background-free" spectra 

The excitation energy of the GMR is expressed in the scaling model [lOj] as: 



e GM r = hJ — Ka (1) 

V m < r z > 

where m is the nucleon mass, < r 2 > is the ground-state mean-square radius, and Ka, 
the incompressibility of the nucleus. In order to determine the incompressibility of infinite 
nuclear matter, K^, from the experimental GMR energies, one builds a class of energy 
functionals, E(p), with different parameters which allow calculations for nuclear matter and 
finite nuclei in the same theoretical framework. The parameter-set for a given class of energy 
functionals is characterized by a specific value of K^. The GMR strength distributions are 
obtained for different energy functionals in a self-consistent RPA calculation. The Koo 
associated with the interaction that best reproduces the GMR energies is, then, considered 



the "correct" value. This procedure, first proposed by Blaizot is now accepted as 



the best way to extract K^ from the GMR data and, following this procedure, it has 
been established that both relativistic and non-relativistic calculations are now in general 
agreement with = 240 ± 10 MeV [l2, 3, Q. 

The determination of the asymmetry term, K T , associated with the neutron excess (N- 
Z), remains very important because this term is crucial in obtaining the radii of neutron 



stars in EOS calculations 
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13, 



181 ] . Indeed, the radius of a neutron star whose 



mass is between about 1 and 1.5 solar masses (Mp) is mostly determined by the density 



dependence of the symmetry-energy term 
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20] . Previous attempts to extract this term 



from experimental GMR data have resulted in widely different values, from -320±180 MeV 
in Ref. (21] to a range of -566±1350 MeV to 139±1617 MeV in Ref. [22]. Measurements 
of the nuclear incompressibility over a series of isotopes provide a way to "experimentally" 
determine this asymmetry term in a direct manner. The Sn isotopes (A=112-124) afford 
such an opportunity since the asymmetry ratio, ((N-Z)/A), changes by more than 80% over 
this mass range. 

In this Letter, we report on new measurements on GMR in the even-A Sn isotopes. 



The GMR has been identified previously in some o: 



measurements with inelastic a-scattering 



flfl 
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the Sn isotopes compact peak in 



241 ] and although resonance parameters 



for GMR in the Sn isotopes close to the values reported here have been extracted in the past 
using less accurate techniques [211 ]. the potentially large systematic errors in those values 
necessitated the present measurements where such problems have been eliminated. We find 
that the GMR energies in the Sn isotopes are lower than the values predicted in recent 
theoretical calculations even though the interactions used in these calculations reproduce 
the GMR energies in the "standard" nuclei, 208 Pb and 90 Zr, very well. Also, we obtain a 
value K T = -550 ± 100 MeV from this data. 

The experiment was performed at the ring cyclotron facility of the Research Center for 
Nuclear Physics (RCNP), Osaka University, using inelastic scattering of 400- MeV a par- 
ticles over the angular range 0°-8.5°. Details of the experimental technique and the data 
analysis procedure have been provided previously 0, o, 8] and are only briefly described 
here. Inelastically-scattered a particles were momentum-analyzed with the high-resolution 



magnetic spectrometer "Grand Raiden" 



251 ] and detected in the focal-plane detector system 



comprised of two multi-wire drift chambers and two scintillators, providing particle identi- 
fication as well as the trajectories of the scattered particles. The vertical position spectrum 
obtained in the double-focused mode of the spectrometer was exploited to eliminate all 
instrumental background ay , B]. The background-free "0°" inelastic spectra for the Sn 
isotopes are presented in Fig. [IJ In all cases, the spectrum is dominated by the GMR peak 
near E x ~ 15 MeV. 

In order to extract the GMR strengths, we have employed the now standard MDA proce- 
dure 26j. The cross-section data were binned into 1-MeV energy intervals between 8.5-31.5 
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FIG. 1: Excitation-energy spectra for all even-A Sn isotopes, obtained from inelastic a scattering 
at 9 lab = 0.69°. 

MeV and for each excitation energy bin, the experimental 17-point angular distribution 
d rffT (^m; E x ) was fitted by means of the least-square method with the linear combination 
of calculated distributions -^-(9 cm , E x ), so that: 

j exp 7 rlrr cal 

^(Oam, E x ) = J2 ME,) x =^-{e m , E x ) (2) 

^^cal 

where -^-(9 cm , E x ) is the calculated distorted-wave Born approximation (DWBA) cross 
section corresponding to 100% energy-weighted sum-sure (EWSR) for the L-th multipole. 
This procedure provides strength distributions simultaneously for various multipoles. 

The DWBA calculations were performed following the method of Satchler and Khoa 



271 ] using density-dependent single folding model, with a Gaussian a-nucleon potential for 



the real part, and a Woods-Saxon imaginary term. We used the transition densities and 



sum rules for various multipolarities as described in Ref. [28f. The optical model (OM) 
parameters were obtained from analysis of elastic scattering cross sections measured in a 
companion experiment. 

Although all strength distributions up to L=3 have been reliably extracted from the 
multipole decomposition, only the GMR strengths, the focus of this paper, are shown in 
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FIG. 2: GMR strength distributions obtained for the Sn isotopes in the present experiment. Error 
bars represent the uncertainty due to the fitting of the angular distributions in MDA. The solid 
lines show Lorentzian fits to the data. 

Fig. [2J The solid lines in the figure represent Lorentzian fits to the observed strength 
distributions. The choice of the Lorentzian shape is arbitrary; the final results are not 
affected in any significant way by using, instead, a Gaussian shape, for example. The finite 
strength at the higher excitation energies is attributable to the mimicking of L=0 angular 
distribution by components of the continuum 4, f|. The extracted GMR-peak parameters 
and the various moment ratios typically used in theoretical calculations are presented in 
TableE 

The moment ratios, mi/m , for the GMR strengths in the Sn isotopes are shown in 



Fig. [3j and compared with recent theoretical results from Colo (non-relativistic) 12|, [29] 
and Piekarewicz (relativistic) 



13l . |3(| . As can be seen, the calculations overestimate the ex- 
perimental GMR energies significantly (by almost 1 MeV in case of the higher-A isotopes). 
This is very surprising since the interactions used in these calculations are those that very 
closely reproduce the GMR centroid energies in 208 Pb and 90 Zr. Admittedly, there are uncer- 
tainties associated with the range over which the experimental and theoretical distributions 
are compared, and also with the assumptions inherent in the calculations regarding widths. 
However, the calculations reported here are identical in all respects to those performed 
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TABLE I: Lorentzian-fit parameters and various moment-ratios for the GMR strength distributions 
in the Sn isotopes, as extracted from MDA in the present work, is the fc-th moment of 
the strength distribution: m& = / E^S(E x )dE x . All moment ratios have been calculated over 
E x = 10.5-20.5 MeV. The errors quoted for EWSR are statistical only. 



Target E GM r (MeV) T (MeV) EWSR m 1 /m Q (MeV) y / m 3 /m 1 (MeV) y Jm 1 /m- 1 (MeV) 



112 Sn 


16.1 ±0.1 


4.0 ±0.4 


0.92 ±0.04 


16.2 ±0.1 


16.7 ±0.2 


16.1 ±0.1 


114 Sn 


15.9 ±0.1 


4.1 ±0.4 


1.04 ±0.06 


16.1 ±0.1 


16.5 ±0.2 


15.9 ±0.1 


116 Sn 


15.8 ±0.1 


4.1 ±0.3 


0.99 ±0.05 


15.8 ±0.1 


16.3 ±0.2 


15.7 ±0.1 


118 Sn 


15.6 ±0.1 


4.3 ±0.4 


0.95 ±0.05 


15.8 ±0.1 


16.3 ±0.1 


15.6 ±0.1 


120 Sn 


15.4 ±0.2 


4.9 ±0.5 


1.08 ±0.07 


15.7 ±0.1 


16.2 ±0.2 


15.5 ±0.1 


122 Sn 


15.0 ±0.2 


4.4 ±0.4 


1.06 ±0.05 


15.4 ±0.1 


15.9 ±0.2 


15.2 ±0.1 


124 Sn 


14.8 ±0.2 


4.5 ±0.5 


1.03 ±0.06 


15.3 ±0.1 


15.8 ±0.1 


15.1 ±0.1 



for 208 Pb and 90 Zr, and the experimental and theoretical centroids reported here have been 
calculated over exactly the same excitation-energy range. This disagreement remains a chal- 
lenge for the theory: Why are the tin isotopes so "soft"? Are there any nuclear structure 
effects that need to be taken into account to describe the GMR energies in the Sn isotopes? 
The incompressibility of a nucleus, Ka, may be expressed as: 

K A ~ K vol (l + cA- l ' z ) + K T ((N - Z)/Af + K Co uiZ 2 A^ 3 (3) 



Here, c ~ — 1 31], and Kcoui is essentially model-independent (in the sense that the de- 
viations from one theoretical model to another are quite small), so that the associated 
term can be calculated for a given isotope. Thus, for a series of isotopes, the difference 
Ka — KcouiZ 2 'A -4 / 3 may be approximated to have a quadratic relationship with the asym- 
metry parameter, of the type y = A + Bx 2 , with K T being the coefficient, B, of the quadratic 



term. It should be noted that it has been established previously [22|, [32J that fits to the 
above equation do not provide good constraints on the value of K^. However, this expres- 
sion is being used here not to obtain a value for K^,, but, rather, only to demonstrate the 
approximately quadratic relationship between Ka and the asymmetry parameter. 

Fig. H] shows the difference Ka — Kc uiZ 2 A~^ 3 for the Sn isotopes investigated in this 
work vs. the asymmetry parameter, ((N — Z) / A) . The values of Ka have been derived using 
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FIG. 3: Systematics of the moment ratios mi/mo for the GMR strength distributions in the Sn 
isotopes. The experimental results (filled squares) are compared with results of non-relativistic RPA 
calculations by Colo [29] (filled circles) and relativistic calculations of Piekarewicz [30( (triangles). 



241 ] are also shown 



Results for 112 Sn, 116 Sn and 124 Sn reported by the Texas A & M group 
(inverse triangles). The differences between the present results and the Texas A & M results for 
H2,i24g n ^ght j-, e attributable to the background subtraction required in their analysis. 



the customary moment ratio ymi/ra_i for energy of the GMR in Eq. (1). A quadratic fit to 
the data is also shown. The fit gives K T = —550 ± 40 MeV, with the uncertainty attributed 
only to the fitting procedure. Including the uncertainties in in the fit adds another ~25 
MeV to this "error" (to ± 67 MeV) and the uncertainty in the value of Kc ou i (± 0.7 MeV; 



see Ref. 



33] ) would contribute ~ 15 MeV. Considering, further, the approximation made 



in arriving at the quadratic expression, the actual total uncertainty would be somewhat 
larger still; hence the rounded value K T = —550 ± 100 MeV quoted earlier in the text. 



This result is consistent with the value K T = —500 ± 50 MeV obtained recently 



analysis of the isotopic transport ratios in medium-energy heavy-ion reactions 



34 



rom an 



shown in Ref. 



35|. As 



181 ]. this value provides constraints on the radius of a 1.4 M & neutron star 



that are in rather good agreement with recent observational data. Thus, from the data 
on the compressional-mode giant resonances, we now have "experimental" values of both 
and K T which, together, can provide a means of selecting the most appropriate of the 
interactions used in EOS calculations. For example, this combination of values for and 
K T essentially rules out a vast majority of the Skyrme-type interactions currently in use in 
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FIG. 4: Systematics of the difference Ka — KcouiZ 2 A 4 / 3 in the Sn isotopes as a function of the 
"asymmetry-parameter" ((N-Z)/A); K Cou i = -5.2 MeV 
quadratic fit to the data. 



33j ] . The solid line represents a least-square 



nuclear structure calculations [33]. A similar conclusion was reached for EOS equations in 



Ref. 36|. Furthermore, a more precise determination of K T provides additional motivation 



for measurement of isoscalar monopo 
at RIKEN and GANIL, for example 



e strength in unstable nuclei, a focus of investigations 



37 



38] 



In summary, we have measured the energies of the isoscalar giant monopole resonance 
(GMR) in the even- A 112 ~ 124 Sn isotopes via inelastic scattering of 400-MeV a particles at 
extremely forward angles, including 0°. The GMR energies are significantly lower than those 
predicted for these isotopes by recent calculations. Further, the asymmetry-term, K T , in the 
expression for the nuclear incompressibility has been determined to be —550 ± 100 MeV. 

We wish to express our gratitude to G. Colo and J. Piekarewicz for providing results of 
their calculations prior to publication. This work has been supported in part by the National 
Science Foundation (Grants No. INT03-42942 and PHY04-57120), and by the Japan Society 
for the Promotion of Science (JSPS). 
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